Towards the simulation of biomolecules: optimisation of peptide-capped glycine using FFLUX

YesThe optimisation of a peptide-capped glycine using the novel force field FFLUX is presented. FFLUX is a force field based on the machine-learning method kriging and the topological energy partitioning method called Interacting Quantum Atoms. FFLUX has a completely different architecture to that of traditional force fields, avoiding (harmonic) potentials for bonded, valence and torsion angles. In this study, FFLUX performs an optimisation on a glycine molecule and successfully recovers the target density-functional theory energy with an error of 0.89 ± 0.03 kJ mol−1. It also recovers the structure of the global minimum with a root-mean-squared deviation of 0.05 Å (excluding hydrogen atoms). We also show that the geometry of the intra-molecular hydrogen bond in glycine is recovered accurately.EPSRC Established Career Fellowship [grant number EP/K005472

Distinguishing epimers through raman optical activity

The Raman optical activity spectra of the epimers β-d-glucose and β-d-galactose, two monosaccharides of biological importance, have been calculated using molecular dynamics combined with a quantum mechanics/molecular mechanics approach. Good agreement between theoretical and experimental spectra is observed for both monosaccharides. Full band assignments have been carried out, which has not previously been possible for carbohydrate epimers. For the regions where the spectral features are opposite in sign, the differences in the vibrational modes have been noted and ascribed to the band sign changes

Calculation of Raman optical activity spectra for vibrational analysis

By looking back on the history of Raman Optical Activity (ROA), the present article shows that the success of this analytical technique was for a long time hindered, paradoxically, by the deep level of detail and wealth of structural information it can provide. Basic principles of the underlying theory are discussed, to illustrate the technique's sensitivity due to its physical origins in the delicate response of molecular vibrations to electromagnetic properties. Following a short review of significant advances in the application of ROA by UK researchers, we dedicate two extensive sections to the technical and theoretical difficulties that were overcome to eventually provide predictive power to computational simulations in terms of ROA spectral calculation. In the last sections, we focus on a new modelling strategy that has been successful in coping with the dramatic impact of solvent effects on ROA analyses. This work emphasises the role of complementarity between experiment and theory for analysing the conformations and dynamics of biomolecules, so providing new perspectives for methodological improvements and molecular modelling development. For the latter, an example of a next-generation force-field for more accurate simulations and analysis of molecular behaviour is presented. By improving the accuracy of computational modelling, the analytical capabilities of ROA spectroscopy will be further developed so generating new insights into the complex behaviour of molecules

The Raman optical activity of β-D-xylose: where experiment and theory meet

Besides its applications in bioenergy and biosynthesis, β-D-xylose is a very simple monosaccharide that exhibits relatively high rigidity. As such, it provides the best basis to study the impact of different solvation shell radii on the computation of its Raman optical activity (ROA) spectrum. Indeed, this chiroptical spectroscopic technique provides exquisite sensitivity to stereochemistry, and benefits much from theoretical support for interpretation. Our simulation approach combines density functional theory (DFT) and molecular dynamics (MD) in order to efficiently account for the crucial hydration effects in the simulation of carbohydrates and their spectroscopic response predictions. Excellent agreement between the simulated spectrum and the experiment was obtained with a solvation radius of 10 Å. Vibrational bands have been resolved from the computed ROA data, and compared with previous results on different monosaccharides in order to identify specific structure–spectrum relationships and to investigate the effect of the solvation environment on the conformational dynamics of small sugars. From the comparison with ROA analytical results, a shortcoming of the classical force field used for the MD simulations has been identified and overcome, again highlighting the complementary role of experiment and theory in the structural characterisation of complex biomolecules. Indeed, due to unphysical puckering, a spurious ring conformation initially led to erroneous conformer ratios, which are used as weights for the averaging of the spectral average, and only by removing this contribution was near perfect comparison between theory and experiment achieved

Consistent patterns of common species across tropical tree communities

Trees structure the Earth’s most biodiverse ecosystem, tropical forests. The vast number of tree species presents a formidable challenge to understanding these forests, including their response to environmental change, as very little is known about most tropical tree species. A focus on the common species may circumvent this challenge. Here we investigate abundance patterns of common tree species using inventory data on 1,003,805 trees with trunk diameters of at least 10 cm across 1,568 locations1,2,3,4,5,6 in closed-canopy, structurally intact old-growth tropical forests in Africa, Amazonia and Southeast Asia. We estimate that 2.2%, 2.2% and 2.3% of species comprise 50% of the tropical trees in these regions, respectively. Extrapolating across all closed-canopy tropical forests, we estimate that just 1,053 species comprise half of Earth’s 800 billion tropical trees with trunk diameters of at least 10 cm. Despite differing biogeographic, climatic and anthropogenic histories7, we find notably consistent patterns of common species and species abundance distributions across the continents. This suggests that fundamental mechanisms of tree community assembly may apply to all tropical forests. Resampling analyses show that the most common species are likely to belong to a manageable list of known species, enabling targeted efforts to understand their ecology. Although they do not detract from the importance of rare species, our results open new opportunities to understand the world’s most diverse forests, including modelling their response to environmental change, by focusing on the common species that constitute the majority of their trees.Publisher PDFPeer reviewe